Nuclear spin dynamics and Zeno effect in quantum dots and defect centers 
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We analyze nuclear spin dynamics in quantum dots and defect centers with a bound electron under 
electron- mediated coupling between nuclear spins due to the hyperfine interaction ( " J-coupling" in 
NMR). Our analysis shows that the Overhauser field generated by the nuclei at the position of 
the electron has short-time dynamics quadratic in time for an initial nuclear spin state without 
transverse coherence. The quadratic short-time behavior allows for an extension of the Overhauser 
field lifetime through a sequence of projective measurements (quantum Zeno effect). We analyze 
the requirements on the repetition rate of measurements and the measurement accuracy to achieve 
such an effect. Further, we calculate the long-time behavior of the Overhauser field for effective 
electron Zeeman splittings larger than the hyperfine coupling strength and find, both in a Dyson 
series expansion and a generalized master equation approach, that for a nuclear spin system with a 
sufficiently smooth polarization the electron-mediated interaction alone leads only to a partial decay 
of the Overhauser field by an amount on the order of the inverse number of nuclear spins interacting 
with the electron. 



I. INTRODUCTION 

Technological advancements have made it pos- 
sible to confine very few electrons in a va- 
riety of nanostructures such as nanowires, 
quantum dots, donor impurities, or defect 
centers .i- 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 ' 13 - 14 ' 15 - 16 ' 17 - 18 - 19 - 20 - 21 - 22 
One driving force behind these achievements is a series 
of proposals for using the spin of an electron as a 
qubit for quantum computingi 23 i 24 i 25 This spin interacts 
with the nuclear spins in the host material via the 
hyperfine interaction. While this interaction leads to 
decoherence of the electron spin state on one hand, it 
also provides the opportunity to create a local effective 
magnetic field (Overhauser field) for the electron by 
inducing polarization in the nuclear spin system, which 
could be used, e.g., for rapid single-spin rotations^. 
Polarizing the nuclear spin system is also one possible 
way to suppress hyperfine-induced decoherence^ 7 - or it 
can be used as a source of spin polarization to generate 
a spin-polarized current. In any case, controlling the 
dynamics of the Overhauser field and, in particular, 
to prevent its decay, is thus of vital importance in the 
context of spintronics and quantum computation^ 

In GaAs quantum dots the Overhauser field can be- 
come as large as 5T. The build-up, decay, and correla- 
tion time of the Overhauser field have been studied in a 
number of systems^^^^^ 5 -^^^^ 4 ^ 2 - sug _ 
gesting timescales for the decay on the order of seconds, 
minutes, or in one case, even hours £2- 

The dynamics of the Overhauser field are governed by 
the mutual interaction between the nuclear spins. There 
is on one hand the direct dipolar coupling between the 
nuclear spins. On the other hand, due to the presence 
of a confined electron, there is also an indirect interac- 
tion: The coupling of the nuclear spins to the electron via 
the hyperfine interaction leads to an effective interaction 



between the nuclear spins that is known as the electron- 
mediated interaction. While the effect of this electron- 
mediated interaction on the decoherence of the electron 
has been studied previously , 44 ' 45 ' 46 theoretical studies of 
the decay of the Overhauser field have so far studied di- 
rect dipole-dipole interaction and the effect of the hy- 
perfine interaction was taken into account through the 
Knight shift that the electron induces via the hyperfine 
interaction^ 7 - In this article we investigate the effect of 
the electron-mediated interaction between nuclear spins 
on the dynamics of the Overhauser field. While the direct 
dipolar coupling is always present, it can be weaker than 
the electron-mediated interaction for magnetic fields that 
are not too large and may be further reduced via NMR 
pulse sequences or by diluting the concentration of nu- 
clear spins4£ We find in our calculation that, for effective 
electron Zeeman splittings u> (sum of Zeeman splittings 
due to the external magnetic field and Overhauser field) 
larger than the hyperfine coupling strength A, the decay 
of the Overhauser field due to the electron-mediated in- 
teraction is incomplete, i.e., that only a small fraction 
of the Overhauser field decays. In a short-time expan- 
sion that is valid for u> larger than A/y/N, where N is 
the number of nuclear spins with which the electron in- 
teracts, we find a quadratic initial decay on a timescale 
T e = N z / 2 lj/A 2 . We show that, by performing repeated 
projective measurements on the Overhauser field, a quan- 
tum Zeno effect occurs, which allows one to preserve the 
Overhauser field even for relatively small effective elec- 
tron Zeeman splittings larger than A/y/N. 

In Sec. |TT] we briefly review the quantum Zeno effect 
and give the corresponding main results for the case of 
the Overhauser field. We start our detailed discussion in 
Sec. IIIII by writing down the Hamiltonian for the hyper- 
fine interaction and by deriving an effective Hamiltonian 
for the electron-mediated interaction. In Sec. IIVI we 
derive an expression for the short-time behavior of the 
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Overhauser field mean value. In Sees. IVl and |VT1 we ad- 
dress the long-time decay of the Overhauser field due to 
the electron-mediated interaction. Some technical details 
are deferred to Appendices E] and [B] 



II. ZENO EFFECT 

The suppression of the decay of a quantum state due 
to frequently repeated measurements is known as the 
quantum Zcno effect. The concept of the quantum 
Zeno effect^ is almost as old as quantum mechanic o 50 i 51 
and it remains one of the most intriguing quantum ef- 
fects. It has been studied intensively from the theoreti- 
cal sid o 52 i 53 ' 54 and also experimental evidence has been 
found in recent years.— 

For a two-level system initialized to the exited state, 
the survival probability P s in the exited state as a func- 
tion of the elapsed time t is initially given by P s (t) = 
1 — c s t 2 /r 2 , with the constant c s and the timescalc t s be- 
ing system dependent. A projective measurement at time 
T m resets the system to the excited state with probability 
Ps(j~ m ). Repeating the measurement m times at inter- 
vals T m <C t s , the survival probability is F s . me ns(mT m ) = 
(l-c s r^/r s 2 ) m « l-c s mr ra /(T s 2 /T ra ), for cmr^/r 2 < 1. 
The survival probability at time t = mT m is thus in- 
creased due to the frequently repeated measurements: 
instead of a quadratic decay on a timescale r s without 
measurements, we have a linear decay on a timescale 
Ts/r m . 

A more complex observable such as the mean of the 
Overhauser field (h z (t)) = Tr{h z p(t)}, may also show 
a Zcno effect. That (h z (t)) shows an initial quadratic 
decay is, however, not obvious and actually depends on 
the initial state of the nuclear spin system pi(0). For 
the short-time behavior of (h z (t)} we expand in a Taylor 
series 

(M*)) = (M0))+*(Mi + ~(Ma + -.., (1) 

with (h z ) n = d n (h z (t))/dt n \ t=() . If (h z ) 1 = 0, the i-linear 
term vanishes and the initial decay is quadratic in time. 
In Sec. IIVI we calculate the initial dynamics of (h z (t)) 
and explain the conditions under which {h z )\ = 0. We 
find an initial decay of the form 
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The timescale r e and the constant c are given below in 
Eq. dTTJ) and Eq. (18]) respectively. 

Let us now consider a sequence of repeated measure- 
ments of the Overhauser field h z (t). In the context 
of quantum dots, several proposal a 56 i 57 ' 58 to implement 
such measurements have been put forward. A mea- 
surement of h z shall be performed after a time r m . If 
this measurement is projective, i.e., if it sets all the off- 
diagonal elements of the density matrix in a basis of h z - 
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FIG. 1: Effect of projective measurements at time intervals 
T m = Te/10 on the time evolution of the Overhauser field 
expectation value (h z (t)). Due to the Zeno effect, the de- 
cay with measurements is 1 — ct/T zeno rather than 1 — ct 2 /t| 
without measurements, where r zmo = T 2 /r m . The formula 
1 — ct/T zeno for the decay with measurement is only strictly 
valid at times t = m% with m being a positive integer. Af- 
ter the measurement at t = mr m the decay is again quadratic 
with time dependence (h z (mrm)) / (h z (0)) — c(t — mr m ) 2 /r,? 
(broken lines). 



cigenstatcs to zero (we discuss requirements on the accu- 
racy of the measurement in Appendix |B|) . the dynamics 
after r m again follow Eq. Repeating the measure- 

ment at times 2T m ,3r m , leads to a change of the 
decay of the Overhauser field in the same way as we de- 
scribed it for the two-level system above: 



{hz(t)h 



(Mo)) 



= l-c- 



t 



— ■ (3) 



Instead of a quadratic decay oc t 2 /r 2 we have a linear 
decay oc t/T zeno with r zeno — t 2 / T m . We note that the 
expression for (h z (t)) zeno in Eq. ((3]) is only strictly valid 
at times mr m with m being a positive integer. Between 
these times (h z (t)} changes according to Eq.®. The 
derivation of Eq. ^ requires citit^/t 2 = ctjr zeno -c 1. 
Fig. [T] shows the Zeno effect, i.e., the difference between 
(M*))/(M0)) and (h z (t)) zeno /(h z (Q)). 

In addition to requirements on the measurement accu- 
racy (see Appendix [BJ), the results in this section rest on 
the following separation of timescales: 



7pm < ^ TV/ 



(4) 



where r pm is the time required to perform a single mea- 
surement and t x the timescale up to which the short- 
time expansion for (h z (t)} is valid. In general, t x can be 
shorter than r e . A specific case (fully polarized nuclear 
state), where the short-time expansion has only a very 
limited range of validity, is discussed in Sec. IIV Al For 
the systems studied in experiment, we expect t x to be 
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comparable to or longer than r e , since the experiments 
performed so far show timescales for the decay of (h z (t)) 
on the order of seconds, minutes, or in one case, even 
hours4£ We note that it may be a demanding task to 
perforin the fast and precise measurements required to 
obtain a Zeno effect in the present content. Still, experi- 
mental progress in the control of the nuclear field, such as 
that shown in Ref. |43|, suggests that such measurements 
may be within reach in the near future. 

We continue our discussion by deriving the effective 
Hamiltonian we use both for calculating short-time dy- 
namics and the long-time behavior of (h z (t)}. 

III. HAMILTONIAN 

We aim to describe the dynamics of many nuclear spins 
surrounding a central confined electron spin in a mate- 
rial with an s-type conduction band (e.g. GaAs, Si, etc.), 
where the dominant type of hyperfine interaction is the 
Fermi contact hyperfine interaction. The electron may 
be confined in many nanostructures such as nanowires, 
quantum dots or defect centers. Under the assumption 
that other possible sources of nuclear spin dynamics, such 
as nuclear quadrupolar coupling, are suppressed^ the 
two strongest interactions between nuclear spins in these 
nanostructures are the electron- mediated interaction ( " J- 
coupling" in NM R 60 i 61 ) and the direct dipolc-dipole in- 
teraction. It turns out that, for a large number of nuclei 
N and up to magnetic fields of a few Tesla (for GaAs), the 
contribution of the electron-mediated interaction to the 
initial decay of the Overhauser field is dominant (see Ap- 
pendix [X]). The Hamiltonian contains three parts: The 
electron and nuclear Zeeman energies and the Fermi con- 
tact hyperfine interaction: 

H = H e + H„ + H en = e z S z +r lz Y J I Z k+S-h. (5) 

k 

Here, the operator 

h = Y,A k T k (6) 

is the Overhauser field. Further, S is the electron spin 
and Jfc the nuclear spin at lattice site k that couples 
with strength A k = AvQ\ij){rk)\ 2 to the electron spin, 
where A = ^ fc A k is the total hyperfine coupling con- 
stant, vq the volume occupied by a single- nucleus unit 
cell and ip(r k ) the electron envelope wave function. We 
define the number of nuclear spins TV interacting with 
the electron as the number of nuclear spins within an 
envelope- function Bohr radius of the confined electron^ 
The Bohr radius o,b for an isotropic electron envelope is 
denned through^ i/)(r k ) = -0(O)e-( rfc / aB )"/ 2 ,where q = 1 
gives a hydrogen-like wave function and q = 2 a Gaus- 
sian. Finally, e z and rj z are the electron and nuclear 
Zeeman splittings, respectively (we consider a homonu- 
clear system). We derive an effective Hamiltonian for 



the electron-mediated interaction between nuclear spins, 
which is valid in a sufficiently large magnetic field. Us- 
ing a standard Schrieffer- Wolff transformation^ H e g = 
e s He~ s , in lowest order in H en , with the transforma- 
tion matrix S = A k ((e z + h z - ij z + A k /2)- 1 S + I^ 
~{e z + h z - r] z - A k /2)~ 1 S-l£) /2, which eliminates 
the off-diagonal terms between electron and nuclear 
spins, we find the effective Hamiltonian iJ c g ~ Hq + V 
(similar to Refs. HHHl), where: 

H a = e z S z +?7 2 ^/£l + S z h Z: (7) 

k 

v = ] ——({h_,h + }s z + hh.,h + ]). 

4(e z - rjz+hz) \ 2 J 

(8) 

In Eq. ([5]) we have neglected terms which are suppressed 
by a factor A k /(e z — rj z + h z ) and the raising and low- 
ering operators are defined as S± = S x ± iS y and sim- 
ilarly for h± and Ij^. The commutator [h—,h+] is de- 
fined in the usual way and {h—, h+} = h-h + + h + h- is 
the anti-commutator of h- and h + . We note that H e g 
neglects the transfer of spin polarization from the elec- 
tron to the nuclei. The electron transfers an amount of 
angular momentum to the nuclear system on the order 
(A/y/Nuj) 2 < I for lo > A/\/N. For to ~ A these contri- 
butions are suppressed by a factor of 0(1 /N) compared 
to the decay of (h z (t)) under H c ff. For very special initial 
states, where H c g leads to no dynamics, e.g., for uniform 
polarization, the transfer of spin from the electron to the 
nuclei is the only source of nontrivial nuclear spin dy- 
namics and therefore should be taken into account. We 
discuss one such initial state, namely, a fully polarized 
nuclear system, in Sec. IIV A[ 

In the following we further replace h z in the denomi- 
nator of Eq. (|8|) by its initial expectation value (h z ) = 
Tr{/i 2 p(0)} and introduce the effective electron Zeeman 
splitting 

oj = e z - Vz + (h z )^e z + (h z ). (9) 

This replacement assumes that the initial state does not 
change significantly and is valid up to corrections sup- 
pressed by ct/uj, compared to the dynamics under H e ff. 
Here a = yf (hi) — (h z ) 2 is the initial width of h z . For 
an unpolarized equilibrium (infinite temperature) nuclear 
spin state we have a cx A/y/N, limiting the range of va- 
lidity to w > A/y/N. Further restricting our treatment 
to I = 1/2 we may write V as 

v = ^ f S > E 4k4#/f + \ E A U S * - p k) ) . 

y k^l k j 

(10) 

where in the sum over k and I the terms k = I are ex- 
cluded. In the next sections we will discuss the dynamics 
of the Overhauser field both at short and at long times 
in the regimes where a perturbative treatment in V is 
appropriate. 
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IV. SHORT-TIME EXPANSION 

With respect to the Zeno effect as discussed in Sec. 
ITT1 our main interest lies in the short-time behavior of 
(h z (t)) (see Eq. ([1])). To calculate (h z )i and (h z ) 2 , we 
expand 

(h z (t)) = Tr{h z exp (-iHt)p(0) exp (iHt)} (11) 

at short times. The first term (h z (0)) = Tv{h z p(0)} 
gives the expectation value at time zero, while the i-linear 
term is proportional to (h z )i = — iTr{h z [H, p(0)]}. Using 
the cyclicity of the trace we find that Tr{h z [H,p(0)]} = 
Tr{[p(0),h z ]H}. Writing p(Q) = p e (Q) <g> pi(0) we have, 
for an initial nuclear spin state pi(0) without transverse 
coherence, [pi(0), h z ] = and thus the t-linear term van- 
ishes. 

To determine the frequency of projective measure- 
ments required to induce a Zeno effect, we are interested 
in (h z ) 2 = -Tr{h z [H, [H,p(0)]}}. We calculate (h z ) 2 be- 
low using the effective Hamiltonian H a g as derived in Sec. 
IIIII The range of validity is limited by higher-order terms 
in the effective Hamiltonian which are proportional to 
(h+h-) n /uj( n+1 \ n = 2,4,... These higher-order terms 
give corrections to (h z }2 which are suppressed by a fac- 
tor (A/y/Nuj) n . Thus the results for (h z (t)) up to 0(t 2 ) 
given below are valid in the regime to ^> A/y/N. Using 
that [h z ,Ho] = 0, we may simplify (h z ) 2 considerably 
and we find for an arbitrary electron spin state: 

(Ma = -■^ s 'Tri{h,\pi(0),h+h-]h+h-}. (12) 

To further simplify, we assume a product initial state of 
the form 

p(0) = p e (0)®pi(0) =p e (0)®kpl k , (13) 

Plk = 1/2 + f k Il- f k = f k (0) = 2(4! (0)). (14) 

For simplicity we restrict our treatment to I = 1/2 and 
thus fk G [-1,1]. With this we find 



1 J2h A l A ^ii h * ® (i+fjIj)Vk-I')h 



ALU 2 



hi 



(15) 

Evaluating the commutators and the trace, we find for 
the decay of the Overhauser field mean value (h z (t)), up 
to corrections of 0(t 4 ), 



(h z (t)} = (M0)> 



' j2 A l A i( A *- A ^~fi)- 

kl 

(16) 

We note that both for uniform coupling constants Ak = 
A/N and for uniform polarization fk = p, Vfc, the i 2 -term 
vanishes. This is, in fact, what one would expect, since 
H c g only leads to a redistribution of polarization and 
both for uniform polarization and uniform coupling con- 
stants, such a redistribution does not affect h z . Rewriting 



the sum in Eq. (|16p we obtain (again up to corrections 
of 0(t 4 )) 



(M*)) 

(Mo)) 



i 



TV 3 / 2 



LO 



A 2 



(17) 



with the numerical factor c only depending on the distri- 
bution of coupling constants through oik — NAk/A and 
the initial polarization distribution fk through 



1 



32iVc 



^ZaWi{ak-m){fk- fi. 



(18) 



hi 



where Co = ^2kf k(Xk - We note that, up to the factor 
c (see Fig. the timescale r e agrees with a previous 
rough estimate^ for the timescale of nuclear-spin dynam- 
ics under the electron-mediated nuclear spin interaction. 
In Table [I] we give r e for a variety of values of the number 
of nuclear spins N and of lo = e z — rj z + (h z ). 
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TABLE I: This Table gives explicit values for the timescale r e 
of the t 2 term in the short-time expansion of (h z (t)). We give 
T e for various values of the number of nuclear spins TV and of 
lo = € z — rj z + (h z ). When lo = A/\/~N we are at the lower 
boundary of tj-values for which the result for r e is valid. The 
parameters used are relevant for a lateral GaAs quantum dot: 
A = 90/ieV, g = -0.4. 

The coupling constants Ak have a different dependence 
on k, depending on the dimension d and the exponent q 
in the electron envelope wave function through^! Ak — 
Aoe~( k / N ^ q/ . For a donor impurity with a hydrogen-like 
exponential wave function we have d = 3, q = 1, d/q = 3, 
whereas for a 2-dimensional quantum dot with a Gaus- 
sian envelope function we have d = 2, q = 2, d/q = 1. 
In Fig. [2] we show the constant c for the case d/q = 1 
and a particular choice of the polarization distribution. 
We give the dependence on d/q in the inset of Fig. [5] 
While c is independent of N for N > 100, it changes 
considerably depending on the initial nuclear spin state, 
which is parameterized by the fk- Since there are nei- 
ther experimental data nor theoretical calculations on 
the shape of the polarization distribution, we assume for 
the curves in Fig. [2] that it has the same shape as the 
distribution of coupling constants Ak, but with a differ- 
ent width, reflected in the number of nuclear spins N p 
that are appreciably polarized. The motivation for this 
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FIG. 2: Numerical prefactor c (given in Eq. (fl8]l) of the t 2 - 
term in the decay of the Overhauser field mean value (h z (t)}. 
While c turns out to be independent (for iV > 100 in the case 
shown according to numerical summation) of the number N 
of nuclear spins within a Bohr radius of the electron envelope 
wave function, it does depend on the type of structure and 
the initial polarization. We show the case of a 2-d quantum 
dot with a Gaussian electron envelope (d/q — 1). The depen- 
dence on the initial polarization is parameterized by N/N p , 
where N p is the number of nuclear spins that is polarized sub- 
stantially (see text). Inset: dependence of c on the ratio d/q 
for N/Np = 1. We see that, e.g., for a donor impurity with a 
hydrogen- like wave function (d/q — 3) the prefactor c is more 
than three orders of magnitude smaller compared to the 2-d 
lateral quantum dot with d/q = 1. 



choice is that if polarization is introduced into the nu- 
clear spin system via electron-nuclear spin flip-flops, the 
probability for these flip-flops is expected to be propor- 
tional to some power of A k /A^. The degree of polariza- 
tion at the center we denote by p € [—1,1]. We may 



thus write f k = pe 



-(k/N p ) 



We see in Fig. [2] that 



c grows monotonically with N/N p , i.e., a localized po- 
larization distribution (N/N p > 1) decays more quickly 
than a wide spread one (N/N p < 1). 

In the context of state narrowing ) 56 i 57 i 58 the short-time 
behavior of the width of the Overhauser field a(t) = 
y/ (h 2 (t)) — (h z (t)) 2 is also of interest. Nuclear spin state 
narrowing, i.e., the reduction of tr, extends the electron 
spin dccoherence time. Repeating the above calculation 
for (h 2 (t)) and using the result for (h z (t)) we find ( up to 
corrections of 0(t 4 )) for the variance of the Overhauser 
field 



a 2 (t) = a 2 (0) [l + c a 



i 1 



(19) 



with the range of validity to > A/\/N, limited by higher- 
order corrections to the effective Hamiltonian as in the 
case of (h z (t)). Here, the dimensionlcss constant c a is 



\ £ ala 2 (a k - ai)(f k - fi)(j k a k + fiai), 



16N C<J0 



kl 



(20) 

where c CT o = J2 k a k(^ ~ fk)- Taking the square-root of 
a 2 (t) and expanding it for c a t 2 /t 2 « 1 we find for the 
width (up to corrections of 0(t 4 )) 



(21) 



<7(i)=<7(0) ( l + C a — 



Thus, also for the width of the Overhauser field the initial 
dynamics is quadratic in time with the same timescale as 
the mean. 



A. Fully polarized case 

In this section we analyze the special case of a fully 
polarized nuclear spin system, where the effective Hamil- 
tonian derived in Sec. IIIII gives no dynamics and thus the 
corrections due to the transfer of polarization from the 
electron to the nuclei become relevant. We thus must 
return to the full Hamiltonian in Eq. (|5|). Using the 
fact that the total spin J z = S z + Y^k 1S a conserved 
quantity, we transform into a rotating frame where the 
Hamiltonian takes the formal 



H' = (i z + h z )S z + i (h+S_ + h_S 



(22) 



with e = e — ?; z . To have any dynamics for a fully 
polarized nuclear spin system (all spins ||)), the ini- 
tial state of the electron must be + with 
Sjj, =/= 0. Since the Iff) part gives no dynamics we consider 
|"0(O)) = |4; TT ■ ■ ■ ?)■ At any later time we may thus 
write 

Hit)) = a(t) wo)} + M*) \t, TT ■ ■ • TUT • • • T) , 

k 

(23) 

with a(0) = 1 and 6^(0) = 0, V7c. The same case was 
studied in Ref. IfJil . However, this study was performed 
from the point of view of electron spin decoherence. For 
the expectation value of (h z (t)}, we find, in terms of a(t) 
and b k (t), 

(h z (t)) = m)\ h z m)) = 4 - £ \b k (t)\ 2 A k , (24) 



where we have used the normalization condition |a(t)| 2 + 
J2 k \bk(t)\ 2 = 1- Using the time-dependent Schroedingcr 
equation idt \ip(t)) = H' \ip(t)), we obtain the differential 
equations for a(t) and b k (t): 

d(t) = l(2e z + A)a(t)- l -J2b k (t)A k , (25) 

k 

hit) = - l -^a(t)- l -(2e z + A-2A k )b k (t). (26) 
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Inserting a power-series Ansatz a(t) = Y^i a t l an d 
bk(t) — J2i &i ^ into these equations and comparing co- 
efficients yieids recursion reiations of the form 



4(1 + 1) 



I y(2e, + a)a«-^ iyE ^ fc) 



- -^ A *^a^-^—(2e z +A-2A k )b« ) 



2(1 + 1) 4(Z + 1) 



(27) 

e 

(28) 



Iterating these recursion relations using that a(0) = 1 
and 6fe(0) = 0,Vfc, we find, neglecting corrections of 
0(0, 



(MO) 

(MO) 



2A ^ fc 



(29) 



For the case of a 2-d quantum dot with Gaussian enve- 
lope wave function, where we have A k = Ae~ k / N /N, we 
find, evaluating A^ by turning it into an integral in 
the continuum limit N 3> 1, (again up to corrections of 
0(t 4 )) 



(MO) 



(30) 



where r c = N/A. To obtain the range of validity for this 
result we go to higher order in t. Again for the case of a 
2-d quantum dot with Gaussian envelope wave function 
we find up to 0(t 4 ), neglecting terms that are suppressed 
by 0(1 /N) in the i 4 -term, 



(MO) 
(MO) 



i / 1 

'-it 



i 

18 



(31) 



Here, T4 = 2viV / y / A(2e z + A). This shows that in some 
cases the higher order terms in the short-time expansion 
can have a considerably shorter timescale. Comparing 
the short-time expansion with a calculation for (S z ) in 
the case of uniform coupling constants 6 ^ suggests that 
the full dynamics contain oscillations with a frequency 
oc e z + A/2, thus limiting the range of validity of the 
short-time expansion to t <C (e z + A/2)^ 1 . 

With this we finish our discussion of the short-time 
dynamics and of the Zeno effect and move on to long- 
time behavior. We first show the results of a Dyson- 
series expansion in Sec. [V] and in Sec. IVII we treat the 
problem using the generalized master equation, showing 
that the Dyson-series expansion gives the leading-order 
contribution in A/ to. 



up to second order in the interaction V. This allows us 
to obtain the full time dynamics of (h z (t)). Since the 
Dyson-series expansion is not a controlled expansion (it 
leads to non-secular divergences in time at higher order) , 
we will only see from the generalized master equation 
calculation in Sec. IVII that the Dyson series result gives 
the correct leading order contribution in A/to. Thus, the 
results in this section are expected to be valid in the 
regime u> A. 

We transform all operators into the interaction picture 
by O = e tH otQ e -iH a t _ j n ^ ne interaction picture we have 
(h z (t)) = Tr{h z p(t)}, with h z = h z since [H ,h z ] = 0. 
Expanding p(t) in a Dyson scries we find 65 

~p(t) = p(0)-i f dt'[V(t'),p(Q)] 



1 dt> f dt"[V(t'),[V(t"),p(0)]}+O(V 3 ), 



where 
V(t) 



Jo 



= p iH ty e -iH t 



(32) 



S z 



g£ e ^-^ 7 + /r . (33 ) 



We assume again the same initial state as in Sec. IIVI and 
thus the term linear in V will drop out under the trace as 
it only contains off-diagonal terms. From the remaining 
two terms we find 



(MO) = (Mo)) 



1 



A\A\(h-h) 



cos 



{A k -Ai)- 



A k - Ar 



- 1 



(34) 



We first verify that this result is consistent with the short- 
time expansion in Sec. IIVI For this we use that A k < 
Aq oc A/N and thus for times t <C t c = N/A we may 
expand the cosine in the above expression, recovering, to 
second order in t, the result in Eq. (fl"6)) . For the full 
time dynamics we note that the sum over cosines leads 
to a decay on a timescale of t c = N/A, since for t > t c 
the different cosines interfere destructively. We illustrate 
this with an example: for a particular choice of the initial 
polarization distribution (d/q = 1 and N p = N) we may 
evaluate the sum in Eq. (|34[) in the continuum limit and 
find 



(MO) 
(MO)) 



1 



8N^ 9{t/Tc) - 



(35) 



The function g(t) is explicitly given by 



V. DYSON-SERIES EXPANSION 

In this section we calculate the expectation value of 
the Ovcrhauscr field (h z (t)) in a Dyson-series expansion 



5(0 



^ i 4 -16i 2 + 64isinQj - 256 sin 2 Q 

' (36) 

with g(0) = and g(t — > 00) = 1. We thus find a power- 
law decay on a timescale r c by an amount of 0(1/N). 
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0.108 




FIG. 3: In this figure we show the ./V-dependence of 1 — 
{hz}stat/{h z (0)}, i.e., the part by which (h z ) decays in units 
of pA 2 /Nu> 2 , in the regime uj 3> A. This plot is for a 3-d de- 
fect center with a hydrogen- like electron envelope (d/q = 3) 
and the initial polarization is parameterized by N/N p = 0.5 
as described in Sec. IIVI For this choice of polarization 
distribution the decay is of 0(1/N). The inset shows the 
full time dynamics of {h z (i))/{h z (0)} as given in Eq. (|35)l for 
d/q = 1, N 3> 1, N/Np = 1. We see that the decay occurs on 
a timescale of r c = N/A. 



Since the sum of cosines in Eq. decays, the remain- 
ing time-independent sum gives the stationary value (up 
to the Poincarc recurrence time 66 ) 



A 



AN 2 



E 

k^tl 



Ctk - Oil 



(37) 



For a system with a large number of nuclear spins N 3> 1 
and a sufficiently smooth polarization distribution, this 
stationary value differs only by a term of 0(1/ N) from 
the initial value, i.e., (h z ) stat / (h z (0)) = 1-0(1/JV). This 
can be seen in Fig. [3J where we show the N dependence 
of l-(h z ) stat /(h z (0)), i- e -: the part by which (h z ) decays. 
The parameters in Fig. [3] are taken for a 3-d defect center 
with a hydrogen-like electron envelope (d/q = 3) and the 
initial polarization N/N p = 0.5 as described in Sec. IIVI 
For this choice of polarization distribution the decay is 
of 0(1 / N). We also find a 0(1/ N) behavior for other 
values of the parameters d/q and N/N p and thus expect 
this to be generally true for a smoothly varying initial 
polarization distribution. The inset of Fig. [3] shows the 
full time dynamics of (h z (t)} as given in Eq. ([35]) for d/q = 
1,7V > 1,N/N P = 1. 

We note that the 4th order of a Dyson series expansion 
gives secular terms (diverging in t) . We thus move on to 
treat the long-time behavior using a master equation ap- 
proach which avoids these secular terms and shows that 
the Dyson series result gives the correct leading-order 
term in A/uj. 



VI. GENERALIZED MASTER EQUATION 

In this section we study the decay of the Overhauser 
field mean value (h z (i)) using the Nakajima-Zwanzig gen- 
eralized master equation (GME) in a Born approxima- 
tion. The results in this section are valid in the regime 
ui 3> A, since higher-order corrections to the Born ap- 
proximation are suppressed by a factor (A/uj) 2 . 

We start from the GME^ 6 . which for P k p(0) = p(0) 
reads 

P k p(t) = -iP k LP k p(t)- f dt'P k Le-' QL ^QLP kP (t'), 
Jo 

(38) 

where L = Lq + Ly is the Liouville superoperator defined 
as (Lq + Ly)(D = [Hq +V,0], The projection superop- 
erator P k must preserve (I k (t)} and we choose it to have 
the form P k = p e (0)Tr e ® P dk (g) /#fe p h (O)Tr/, where P dk 
projects onto the diagonal in the subspace of nuclear spin 

k and is defined as P dk O = S s =t,i l s ( s fel ( s *l G \ s k)- 
Further, Q = 1 — P k . In a standard Born approximation 
and using the same initial conditions as above, i.e., a 
product state and no transverse coherence in the nuclear 
spin system, we obtain the following integro-differential 
equation for (I k (t)) 



m)) 



3. 

8lu 2 



dT 



cos 



;(A k -Ai) 



r)) - {If (0)}). (39) 



The Born approximation goes to order Ly in the expan- 
sion of the self-energy. Higher-order corrections in Ly are 
estimated to give contributions to the right-hand side of 
Eq. ([39]) that are suppressed by a factor (A/lo) 2 . We 
expect the results of this section to be valid at least for 
lu ^> A, although it could in principle happen that (as 
in the case of the decay of (S z (t))2L) the result for the 
stationary value has a larger regime of validity. On the 
other hand it can not be generally excluded that higher- 
order contributions could dominate at sufficiently long 
times. Integrating Eq. (|39[) we find the formal solution 



MEM) 



o 

z /V 



-(A k ^A l )\ ((P k (t'-r))-(I[(0))). 

(40) 

This shows that (If (i)) = (If (0)) + 0((A/uj) 2 ) and we 
may thus iterate this equation and replace (I k (f — r)) in 
the integral by (i|(0)). This implies, up to corrections 

fO((AM 4 ), 



{im) = {im) 



A 



16uj 2 

l,ljik 

t pt' 

x I dt' drcos 
o Jo 



J2 A f(fk-fi) 



-(A k -A{) . (41) 
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Performing the integrals and summing over the (II (t)) 
weighted by their coupling constants A kl we recover the 
Dyson series result in Eq. (j34|) . This shows that the 
Dyson series expansion gives the leading-order contribu- 
tion in A/lu. 

For the analytical solution of Eq. (|39|) in the station- 
ary limit we perform a Laplace transformation, solve the 
resulting equation in Laplace space, and calculate the 
residue of the pole at s = which yields (up to the re- 
currence time) 

<I fe *) 8tat = lim i f{I z k {t))dt = lhr , s {P k (s)) 

T— too 1 Jq s— »0 

= ^E^(0< J f(* = )). ( 42 ) 

with Z k = We see that (I k ) s tnt is determined 

by weighting the neighboring (If(t = 0)) with the prob- 
ability distribution P^{f)/Zk, which is explicitly given by 

f Aj/iAu-Aif : l^k, 
P k (l) = { (43) 
[ 2cu 2 /A 2 k : l = k. 

We point out that (/|) s tat can be either smaller or larger 
than (If (t = 0)) and that £ fc (^) st at - EfcWC* = «)) 
since the total spin is a conserved quantity. Again ex- 
panding the result in Eq. (|4"2"|) to leading order in A/co 
and summing over the nuclear spins weighted by their 
coupling constants A kl we recover the same result found 
in the Dyson series calculation in Eq. (|37|) . Intuitively 
one would expect a decay even at high fields (although a 
very slow one) to a state with uniform polarization. The 
fact that our calculation shows no such decay suggests 
that the Knight-field gradient, i.e., the gradient in the 
additional effective magnetic field seen by the nuclei, due 
to the presence of the electron, is strong enough to sup- 
press such a decay if the flip-flop terms are sufficiently 
suppressed. 

As discussed in Sec. [V] the decay to the stationary 
value occurs on a timescale r c = N/A. Performing a pro- 
jective measurement at a time t > r m resets the initial 
condition and thus again a small decay occurs. Repeat- 
ing these measurements at intervals longer than r m thus 
allows for a decay of (h z (t)) to zero. 



that is prolonged by a factor r e /r m , where r m is the 
time between consecutive measurements. In Sees. [Viand 
rVTl we have addressed the long-time decay of (h z (t)) us- 
ing a Dyson series expansion and a generalized master 
equation approach. Both show that (h z (t)) only decays 
by a fraction of 0(1/N) for a sufficiently smooth po- 
larization distribution and large magnetic field. It re- 
mains a subject of further study beyond the scope of this 
work whether, and on what timescale, the combination 
of electron-mediated interaction and direct dipole-dipolc 
interaction may lead to a full decay of the Overhauser 
field. Another interesting question concerns the distri- 
bution of nuclear polarization within a quantum dot or 
defect center and its dependence on the method that is 
used to polarize the system. 

We thank G. Burkard, A. Imamoglu, T. Mcunicr, K. 
C. Nowack, D. Stepanenko, J.M. Taylor, M. Trif, and 
in particular F.H.L. Koppens and L.M.K. Vandersypen 
for useful discussions. We acknowledge financial support 
from JST ICORP, the NCCR Nanoscience and the Swiss 
NSF. 



APPENDIX A: ESTIMATION OF 
DIPOLE-DIPOLE CONTRIBUTION 

In this appendix we estimate the timescale arising 
from the direct secular (terms conserving I z .tot = Sfc ^k) 
dipolc-dipole interaction in the short-time expansion of 
the Overhauser field mean value (h z (t)). This gives us the 
range of validity of our calculation in the main text that 
only took into account the electron-mediated interaction 
between nuclei. Let us thus consider the situation where 
the external magnetic field is very high, such that the 
electron-mediated flip-flop terms are fully suppressed. In 
this case the Hamiltonian has the form Hdd = Ho.dd+Vdd, 
with 

k kjtl 

(Al) 

v dd = Ys^tK- ( A2 ) 



VII. CONCLUSION 

We have studied the dynamics of the Overhauser field 
generated by the nuclear spins surrounding a bound 
electron. We focused our analysis on the effect of the 
electron-mediated interaction between nuclei due to the 
hyperfine interaction. At short times we find a quadratic 
initial decay of the Overhauser field mean value (h z (t)} on 
a timescale r e = N z / 2 uj / A 2 . Performing repeated strong 
measurements on h z leads to a Zcno effect with the de- 
cay changing from quadratic to linear, with a timescale 



Here, b k i = ~fj (3 cos 2 (6 k i) - l)/r%, with 9 k i being the 
angle between a vector from nucleus k to nucleus I and 
the z-axis and r k i being the distance between the two 
nuclei^ Further, 7/ is the nuclear gyromagnetic ratio. 
For the short-time expansion, only the off-diagonal terms 
are relevant, since [h z ,H ] = [p(0),Ho] = 0. These off- 
diagonal terms in the case of the electron-mediated in- 
teraction are S z J2k^i A k AiI k If /2tu (see Eq. (TTUjQ . Re- 
placing A k Ai/2cu by b k i in the result for the short-time 
expansion in Eq. (|16|) and also taking into account the fac- 
tor of 1/4 that comes from S 2 in the electron- mediated 
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case wc find 



(ft*(*))dip- 



dip 



kl 

(A3) 

To estimate, we restrict the sum to nearest neighbors 
as the bki fall off with the third power of the distance 
between the two nuclei. Assuming f k = (A k /A ) N / Np wc 
find up to corrections of 0{t A ) 



(h z (t)) c \ip- 



d i p 



<M0)> 



1 - 



i 1 



Td = 



(A4) 



with b being the nearest-neighbor dipole-dipole coupling. 
For GaAs we have b ~ 10 2 s _1 (with 7/ w 10 MHz/TSI). 
For NN P > 1 we have Td ^> 10 _2 s. In the magnetic field 



range shown in Table [I] we thus have S> T e /s/c, which 
justifies neglecting the direct dipole-dipole coupling in 
the short-time expansion. 



APPENDIX B: MEASUREMENT ACCURACY 

The description of the Zeno effect in Sec. [II] relied 
on the assumption that the measurements on h z set 
all off-diagonal elements of the density matrix to zero. 
This assumption requires on one hand a perfect mea- 
surement accuracy for h z (we discuss deviations from 
that below), but on the other hand it also requires the 
/i 2 -eigenstates to be non-degenerate. For non-degenerate 
h z eigenstates a measurement of h z fully determines the 
polarization distribution and we may thus write pi 
after the measurement again as a direct product with 
Pi k { T m) = 1/2 + fk(i~ m )Ik- After the measurement, wc 
thus again have the same time evolution for (h z (t)) as 
given in Eq. (fl6)l . but with ft replaced by /fe(r TO ). Iterat- 
ing Eq. ([TBI) for the case of m consecutive measurements 
at intervals T m one obtains Eq. ([3]). 

Instead of the idealized assumption of a projective 
measurement we now allow for imperfect measurements. 
To describe these measurements we use a so-called 
POVM (positive operator valued measure) £i In a general 
POVM measurement the density matrix changes accord- 
ing to«i 



P 



FyPVFydy, 



(Bl) 



For an unpolarizcd equilibrium (infinite temperature) 
state, the width is a oc A/y/N. Here, a can take any 
value. Let us now consider a measurement that deter- 
mines the value of h z up to an accuracy r\ with a Gaussian 
lincshape. We refer to r\ as the measurement accuracy. If 
the outcome of the measurement is (h z )+y, the diagonal 
of the nuclear spin density matrix after a measurement 
has the form 



Pi{n,n;y) 



2717] 



cxp 



(ft™ - (h z ) - yf 



2if 



(B3) 



Since we aim to describe measurements that at least par- 
tially project the nuclear spin state, we have rj < a. The 
POVM that describes such a measurement is given by 



(B4) 



with 

f(n,y) 



: CXp 



{K - (h z ) - y) 



~ 1 



r?VM^ - V 2 ) 

(K - (K)) 2 



x exp 



2r] 2 
,.2 



2a 2 



2(a 2 -rf)_ 



(B5) 



We note that for 77 <C a we have f(n,y) rts exp(— (ft™ — 
(h z )— y) 2 /2?7 2 )/v / 27rf?. With /(n, y), the operators F y are 
fully determined and it is straightforward to calculate the 
probability for obtaining the measurement result (h z ) +y 



P(y) = 



V /27T((7 2 - if 



■. exp 



2(f7 2 - f] 2 ) 



(B6) 



Clearly, the probabilities add up to one (J P(y)dy = 
1) as they should. Also, when weighting the 
p'j(n,n;y) with their probabilities for occurring, we find 
/ p'j(n,n;y)P(y)dy = pi(n,n). Using Eq. (|B1|) we thus 
find for the matrix elements after a measurement, when 
averaging over all possible measurement outcomes 



Pi(n,m) -> p'j(n,m) = pi(n,m) J \J f(n,y)f(m,y)dy, 

(B7) 

with (for 77 <C er) 



f{n, y) 



1 



7T7) 



■ exp 



when averaging over all possible measurement outcomes 
y. The probability to measure outcome y is given by 
P(y) = Tr{pF y } and the condition J dyF y = 1 ensures Again, for r\ <C c, we thus have 
that the probabilities sum to unity. We consider the nu- 
clear density matrix pi in a basis of h z eigenstates |n) 
with h z \n) = ft™ \n). Wc denote the matrix elements of 
pi by pi(n, m) = (n\ pi \m). For the following description 
we assume that the diagonal of the nuclear spin density 
matrix before the measurement is Gaussian distributed 
around its mean value (h z ) with a width cr, i.e., 



(ft" - h z0 - y) 



21 



2r] 2 



p'i(n,m) = pi(n,m) exp 



(ft™ - ft™) 



m\2 " 



8?y 2 



(B8) 



(B9) 



pi(n,n) 



1 



2-KG 



exp 



{K-{h z )f 

2a 2 



(B2) 



To reduce the off-diagonal elements, the measurement ac- 
curacy must be better than the difference in eigenvalues. 
In the limit r\ — > a projective measurement is recov- 
ered, which sets all off-diagonal elements to zero. Up to 
t 2 in the short-time expansion, only off-diagonal elements 
between states that differ at most by two flip-flops can 
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become non-zero. Thus, to have at least a partial Zeno 
effect^ resulting from the off-diagonal elements being 
partially reduced, the requirement on the measurement 
accuracy is 77 < ft" — h™ with \n) = \m). 
For coupling constants = Ae~ k l N /N , we have typ- 
ically ft" — ft™ oc A/N 3 / 2 . Besides destroying the off- 



diagonal elements of pi through a measurement, there are 
also "natural" dcphasing mechanisms, such as inhomogc- 
neous quadrupolar splittings, electron-phonon coupling, 
or spin-lattice relaxation, that can lead to a reduction of 
the off-diagonal elements of pi . 
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